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Abstract 



Wc consider the linear growth and fragmentation equation 

d d f°° 

g^u(x, t) + — (t(x)u) + j3(x)u = 2 J P(y)K(x, y) u(y, t) dy, 

with general coefficients r, j3 and n. Under suitable conditions (see |16]L the first eigenvalue repre- 
sents the asymptotic growth rate of solutions, also called fitness or Malthus coefficient in population 
dynamics ; it is of crucial importance to understand the long-time behaviour of the population. 
We investigate the dependency of the dominant eigenvalue and the corresponding eigenvector on 
the transport and fragmentation coefficients. We show how it behaves asymptotically as transport 
dominates fragmentation or vice versa. For this purpose we perform suitable blow-up analysis of the 
eigenvalue problem in the limit of small/large growth coefficient (resp. fragmentation coefficient). 
We exhibit possible non-monotonic dependency on the parameters, conversely to what would have 
been conjectured on the basis of some simple cases. 
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Introduction 



Growth and division of a population of individuals structured by a quantity supposed to be conserved 
in the division process may be described by the following growth and fragmentation equation 

d d f°° 

—u{x,t) + —(T(x)u(x,t))+(3(x)u(x,t) = 2j f3(y)K(x,y)u(y,t)dy, x ^ 0, 

u(x,0) = u (x), 
u(0,t) = 0. 

This equation is used in many different areas to model a wide range of phenomena: the quantity 
u(t,x) may represent a density of dusts ([17]). polymers [8j[9], bacteria or cells [HE]; the structuring 
variable x may be size ( |31] and references), label ([2[ [JJ), protein content (P2J [26]), proliferating 
parasite content ([3]); etc. In the literature, it is refered to as the "size-structured equation", "growth- 
fragmentation equation', "cell division equation", "fragmentation-drift equation" or yet Sinko-Streifer 
model. 

The growth speed r = ^ represents the natural growth of the variable x, for instance by nutrient 
uptake or by polymerization, and the rate (3 is called the fragmentation or division rate. Notice that if 
r is such that - is integrable at x = 0, then the boundary condition u(0, t) = is useless. The so-called 
fragmentation kernel n(x, y) represents the proportion of individuals of size x < y born from a given 
dividing individual of size y; more rigorously we should write n(dx,y), with n(dx,y) a probability 
measure with respect to x. The factor "2" in front of the integral term highlights the fact that we 
consider here binary fragmentation, namely that the fragmentation process breaks an individual into 
two smaller ones. The method we use in this paper can be extended to more general cases where the 
mean number of fragments is no > 1 (see [16]). 

Well-posedness of this problem as well as the existence of eigenelements has been proved in [18\ 116], 
Here we focus on the first eigenvalue A associated to the eigenvector U defined by 

d f°° 

— (t(x)U(x)) + {(3{x) + X)U{x) = 2 / p(y) K {x, y)U(y)dy, x > 0, 
& x Jx ^2) 

T U(x = 0) = 0, U{x) > for x > 0, / °° U{x)dx = 1. 

The first eigenvalue A is the asymptotic exponential growth rate of a solution to Problem (TjQ) (see 
It is often called the Malthus parameter or yet the fitness of the population. Hence it 
is of deep interest to know how it depends on the coefficients: for given parameters, is it favorable 
or unfavorable to increase fragmentation ? Is it more efficient to modify the transport rate r or to 
modify the fragmentation rate /3 ? Such concerns may have deep impact on therapeutic strategy 
(see [U [5j [TU [10]) or on experimental design of devices such as PMC A (see [25] and references 
therein). Moreover, when modeling polymerization processes, Equation (TJJ) is coupled with the density 
of monomers V(t), which appears as a multiplier for the polymerization rate (i.e., t(x) is replaced 
by V(t)T(x), and V(t) is governed by one or more ODE - see for instance [22j [231 IS])- Asymptotic 
study of such polymerization processes thus closely depends on such a dependency (see [3 [8] , where 



1 PMCA, Protein Misfolded Cyclic Amplification, is a protocole designed in order to amplify the quantity of prion 
protein aggregates thanks to periodic sonication pulses. In this application, u represents protein aggregates density and 
x their size; the division rate /3 is modulated by sound waves. See Section [3.31 for more details. 
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asymptotic results are obtained under the assumption of a monotonia dependency of A with respect 
to the polymerization rate r). 

Based on simple cases already studied (see [22], [23J [5U US] ) , the first intuition would be that the 
fitness always increases when polymerization or fragmentation increases. Nevertheless, a closer look 
reveals that it is not true. 

To study the dependency of the eigenproblem on its parameters, we depart from given coefficients 
r and f3, and study how the problem is modified under the action of a multiplier of either the growth 
or the fragmentation rate. We thus consider the two following problems: firstly, 



a— (T(x)U a (x)) + {Pix) + \ a )U a {x) = 2 



P(y)K{x,y)U a (y)dy, 



x^O, 



(3) 



rU a (x = 0) = 0, U a {x) > for x > 0, f^U a (x)dx = 1, 



where a > measures the strength of the polymerization (transport) term, as in the Prion problem 
(see [22]), and secondly 



( d f°° 

— ( T (x)V a (x)) + (a/3(x) + A a )V a (x) = 2a J p(y) K (x, y)V a (y)dy, x ^ 0, 

T V a (x = 0) = 0, V a (x) > for x > 0, / °° V a (x)dx = 1, 



(4) 



where a > modulates the fragmentation intensity as for PMCA or therapeutics applied to the cell 
division cycle (see discussion in Section [3]) . 



To make things clearer, let us give some enlightments on the dependency of A a and A a on their 
respective multipliers o and a. First of all, one suspects that if a vanishes or if a explodes, since 
transport dominates, the respective eigenvectors ti a and V a tend to dilute, and on the contrary if a 
explodes or if a vanishes, since fragmentation dominates, they tend to a Dirac mass at zero (see Figure 
[T]for an illustration). But what happens to the eigenvalues X a and A ? Integrating Equation (jl]), we 
obtain the relation 

/•oo 

A a = a /3(x)V a (x) dx 
Jo 

which could give the intuition that A a is an increasing function of a, what is indeed true if (3(x) = (3 is 
a constant since we obtain in this case A a = f3a. However, when /? is not a constant, the dependency 
of the distribution V a (x) on o comes into account and we cannot conclude so easily. A better intuition 
is given by integration of Equation against the weight x, that gives 



A, 



and as a consequence we have that 



J xV a (x) 



dx = r(x)V a (x)dx 



„ t(x) , t(x) 
inf < A a < sup-^. 



x>0 X 



x>0 X 



This last relation highlights the link between the first eigenvalue A a and the growth rate t(x), or more 
precisely l£Ei. For instance, if is bounded, then A a is bounded too, independently of a. Notice 
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Figure 1: Eigenvectors U a (x) for different values of a when f3(x) = 1, k(x, y) = ^o<x<y and r(x) = x. 

In this case we have an explicit expression U a (x) = 2y / a ^/ax + ^f-^ exp (^fax + s §-^. One sees 
that if a vanishes, U a tends to a Dirac mass, whereas it dilutes when a — > +00. 

that in the constant case f3{x) = f3, there cannot exist a solution to the eigenvalue problem ([2]) for 
bounded since we have A a = Ba which contradicts the boundedness of -. In fact we check that 

X 11 ^ X 

existence condition (|26|) in Section 11.21 imposes, for (3 constant, that - is integrable at x = and so 
cannot be bounded. 

X 

Similarly, concerning Equation ([3]), an integration against the weight x gives 



that could lead to the (false) intuition that X a decreases with a - what is indeed true in the limiting 
case t{x) = x. A simple integration gives more insight: it leads to 



This relation makes appear the fragmentation rate j3 when the parameter a is in front of the transport 
term. Moreover, we have seen that when the growth parameter a tends to zero for instance, the 
distribution U a {x) is expected to concentrate into a Dirac mass in x = 0, so the identity X a = 
J (3(x)U a (x) dx indicates that A Q should tend to /3(0). Similarly, when a tends to infinity, X a should 
behave as /3(+oo). 

These intuitions on the link between - and A a on the one hand, f3 and X a on the other hand are 
expressed in a rigorous way by the following result. The main assumption is that the coefficients t(x) 
and (3(x) have power- like behaviours in the neighbourhood of L = or L = +00, namely that 



1 J TU a dx 



a J xU a dx' 



X, 




inf (3(x) < X a < sup(3(x). 

x>0 x>0 



3 76K such that t(x) 



~ TX 

x^tL 



v 



(3(x) 



x^tL 



fix" 1 . 



(5) 



4 



Theorem 1. Under Assumption ([5]), Assumptions (llip -([T2l) on k, and Assumptions (I2ip -(|27 j) given 
in fWjl to ensure the existence and uniqueness of solutions to the eigenproblems ([3]) and we have, 
for L = or L = +oo, 

t( x ) 

lim \ a = lim and lim A a = lim . 

a-¥L x^tL a— >L x ~^t X 

This is an immediate consequence of our main result, stated in Theorem[2]of Section [TTT1 As expected 
by the previous relations, for Problem ([3]) the eigenvalue behavior is given by a comparison between 
f3 and 1 in the neighbourhood of zero if polymerization vanishes (a — > 0), and in the neighbourhood 
of infinity if polymerization explodes (a — > oo). For Problem (j4]), it is given by a comparison between 
r and x (in the neighbourhood of zero when a — > oo or in the neighbourhood of infinity when a — > 0). 

One notices that these behaviors are somewhat symetrical: indeed, the first step of our proof is to 
use a properly-chosen rescaling, so that both problems ([3]) and ([J]) can be reduced to a single one, 
given by Equation (|17p . Theorem [2] studies the asymptotic behaviour of this new problem, what allows 
us to quantify precisely the rates of convergence of the eigenvectors toward self-similar profiles. 

A consequence of these results is the possible non-monotonicity of the first eigenvalue as a function 
of a or o. Indeed, if lim a; _ j .o P( x ) = li m x->-oo P(x) = 0, then the function a *— > X a satisfies lim a _>.o = 
linic^oo X a = and is positive on (0, +oo) because A Q = j f3U a > for a > 0. If lim^o = 
lim^^oo = 0, we have the same conclusion for a i— )■ A a (see Figure [2] for examples). 




1 2 3 4 5 6 7 8 012345678910 



(a) a^Aa (b) a i-> A a 



Figure 2: The dependencies of the first eigenvalue on polymerization and fragmentation parameters 
for coefficients which satisfy the assumptions of Theorem Q] are plotted. The coefficients are chosen to 



obtain non monotonic functions. |2(a 



lima- 
and n(x,y) 



lim T _ 



,/3(x) 



t(x) = 
0, so liniQ^o 



ax"-* 



15+x' 



jr and k(x, y) = ^lo<x<y We have 



1.2x 



y%<x<y- 



We have lim 



x->0 



t(x) 



= lim^oo X a = 0. [2(b)) r(x) - 1+2x $a 
lim^oo ^ = 0, so lim n ^o A a = lim ->-oo A n 



4i 2 
10+x°- s 

0. 
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This article is organised as follows. Section Q] is devoted to state and prove the main result, given 
in Theorem [21 We first detail the self-similar change of variables that leads to the reformulation of 
Problems ([3]) and ([!]) in Problem (]17p . as stated in Lemma [TJ We then recall the assumptions for the 
existence and uniqueness result of [16]. We need these assumptions here not only to have well-posed 
problems but also because the main tool to prove Theorem [2] is given by similar estimates than the 
ones used in [16] to prove well-posedness. 

In Section [21 we give more precise results in the limiting cases, i.e. when lim^^x or lim x ^,L ~^r- 
is finite and positive, and conversely more general results under assumptions weaker than Assump- 
tion ©. 

Finally, Section [3] proposes possible use and interpretation of the results in various fields of applica- 
tion. 



1 Self-similarity Result 

1.1 Main theorem 

The main theorem is a self-similar result, in the spirit of [18]. It declines four times, as Equation ([3]) 
or (UD is considered, and as parameter a or a goes to zero or to infinity. It gathers the asymptotics 
of the eigenvalue and possible self-similar behaviors of the eigenvector, when r and (3 have power- 
like behavior in the neighbourhood of or +oo. We first explain in full details how the study of both 
Equations ([3|) and Q reduce to the study of the asymptotic behavior of a unique problem, as LemmaQ] 
states. 

When fragmentation vanishes or polymerization explodes, one expects the eigenvectors lA a and V a 
to disperse more and more. When fragmentation explodes or polymerization vanishes on the contrary, 
we expect them to concentrate at zero. This leads to the idea of performing an appropriate scaling of 
the eigenvector (U a or V a ), so that the rescaled problem converges toward a steady profile instead of 
a Dirac mass or a more and more spread measure distribution. 

For given k and we define v a or w a by the dilation 

v a (x) = a k U a (a k x), w a (x) = a l V a (a l x). (6) 
The function v a satisfies the following equation 

oo 

( r ( afcx K( x )) + ( A « + Pi ?*)) v *( x ) = 2 J l3{a k y)K{a k x, a k y)v a {y)dy, (7) 

X 

and similarly the function w a satisfies 

oo 

a ~ l ~dx~ (y^)^)) + ( A " + fl / 3 ( a ' x )) w *( x ) = 2o / /3(a l y)n(a l x, a l y)w a (y)dy. (8) 

X 

A vanishing fragmentation or an increasing polymerization will lead the mass to spread more and 
more, and thus to consider the behavior of the coefficients r, (3 around infinity; on the contrary, a 
vanishing polymerization or an infinite fragmentation will lead the mass to concentrate near zero, and 
we then consider the behavior of the coefficients around zero. This consideration drives our main 
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assumption ([5]) on the power-like behavior of the coefficients t(x) and (3{x) in the neighbourhood of 
L = and L = +00. We remind this assumption here 

3 v, 7 G R such that r(x) ~ n", /3(x) ~ /3x 7 . © 

2— »L x^tL 

In the class of coefficients satisfying Assumption some of the assumptions which appear in Sec- 
tion 11.21 to ensure the existence of eigenelements can be written in terms of the power coefficients 7 
and v. Under Assumption ([5]) at L = 0, Assumption (126j) . which is a condition at x = 0, is verified if 
and only if 

7 + 1 - v > 0, (9) 
and in Assumption (|25p in the same way, the condition linking 70 to t(x) becomes 

for v corresponding to L = in Assumption ([5]), 70 + 1 — v > 0. (10) 

In the same way, under Assumption ([5]) at L = +00, Assumption ([271) is also equivalent to Assumption 
([9]). Hence, under both assumptions (f26|) and ([27|) . the coefficients 7 and 1/ in Assumption ([5]) are such 
that 1 + 7 — 1/ > either for L = or for L = +00. 
To preserve the fact that k is a probability measure, we define 

K a (x,y) := a k K(a k x,a k y), K a (x,y) := a l K(a!x,a!y). (11) 

For our equations to converge, we also make the following assumption concerning the fragmentation 
kernels K a and K a : 

For k > 0, 3 k-l s.t. \/ip G C£°(R + ), / ip{x)n a (x,y) dx )• / (p(x)KL(x,y) dx a.e. (12) 

It is nothing but the convergence in a distribution sense of n a (.,y) for almost every y. This is true for 
instance for fragmentation kernels which can be written in an homogeneous form as n(x,y) = -«(|). 
In this case n a is equal to k for all a, so kl = n suits. 

Under Assumption (0), in order to obtain steady profiles, we define 

r a {x) := a~ ku T(a k x), r B (x) := o _,,/ t(o'x), /3 a := a" fc7 /3(a fc x), /3 n := a~ h '/3(a l 'x), (13) 

so that, if A; > 0, there are local uniform convergences on ]R^_ of cind. j3 a : T a — > tx u and 

a— >L 

f3 a — > fix 1 . Equations ([7J) and ([8]) divided respectively by a fc7 and a'^ -1 ^ can be written as 

a->L 

00 

d 



d_ 

dx 



a^-^ — fr^ix)) + (u- k U a + fJ a ) v a (x) = 2 / p a (y)Ka(x,y)v a (y)dy, (14) 



00 



(r a w a (x)) + (a'C 1 -"^ + o 1 -^"- 1 ^^) Mx) = 2a 1 - 1 ^- 1 -^ / /3 a (y)K a (x, y)w a (y)dy. (15) 



X 



In order to cancel the multipliers of r Q and f3 a , it is natural to define 



k = > 0, l = -k = < 0, (16) 

1 + 7 — v 1 + 7 — v 
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which leads to 



d_ 

dx 



{T a (x)v a (x)) + (0 a + p a (x))v a (x) = 2 J P a (y)K a (x,y)v a (y)dy, (17) 



and 



d 

(r a (x)w a (x)) + (G a + P a (x))w a (x) = 2 I (3 a (y)K a (x, y)w a (y)dy, (18) 



dx 



with 



— 7 v — 1 

a^-^c, 6 a = a 1+ T-" A a . 



The signs of k and I express the fact that for a > 1 or a < 1, v a and w a are contractions of t/ aja , 
whereas for a < 1 or a > 1 they are dilations - what is in accordance with our intuition of the 
respective roles of polymerization and fragmentation. Moreover, one notices that if we define a := — , 
then a 1 = a k , so Equations (|17|) and f)18|) are identical. By uniqueness of a solution to this eigenvalue 
problem, it implies that 9 a = Bj_ and v a = wi. We are ready to state this result in the following 

a. a 

lemma. 

Lemma 1. Eigenproblems §3§ and (JU) are equivalent to the eigenproblem (|17j) with k defined by (|16|) . 

t q defined by (|13p . k q defined by (jlip , a = — and i/ie following relations linking the different 
problems: 

v a (x) = a h U a (a k x) = a k Vi(a k x), 9 a = a 1 ^^ 7 X a = a 1 ^^ A i . (19) 

a: a 

Defining (voo,#oo) as the unique solution of the following problem 

( d f°° 

— (tx v v 00 (x)) j t(I3x 1 + 9 00 )v 00 {x) = 2J Py 1 K L (x,y)v oc (y)dy, x ^ 0, 

TVoo(x = 0) = 0, Voo(x) > for x > 0, / °° Voo(x)dx = 1, #00 > 0, 



(20) 



we expect a to converge towards 9^ > and v a towards Voo when a tends to L, so the expressions 
of A Q , A given by (119p will provide immediately their asymptotic behavior. It is expressed in the 
following theorem. 

Theorem 2. Let r, f3 and k verifying Assumptions (|2i p -(|27 p . Let L = or L = +00, and r and 
P verifying also Assumption §5$) (so condition ([9]) is verified). Let n a defined by (jlip with k defined 
by (|16p verifying Assumption (|12|) . Let (t» Q >#«) 6e t/ie unique solution to the eigenproblem (|18p . We 
/iai>e £/ie following asymptotic behaviors 

x r v a (x) >■ x r Voo(x) strongly in L 1 for all r > 0, and 9 a > 9^. 

ot^rL a—tL 

Theorem Q] announced in introduction follows immediately from Theorem [2] and the expression of 
A Q and A a given by (fT9l) . 



S 



1.2 Recall of existence results ([271 OH) 

We first recall the assumptions of the existence and uniqueness theorem for the eigenequation ([2]) (see 
|16j for a complete motivation of these assumptions). This also ensures well-posedness of Problems 
© and g]). 

For all y > 0, n(.,y) is a nonnegative measure with a support included in [0,y]. We define At on 
(1R+) 2 as follows: n(x,y) = for x > y. We assume that for all continuous function ip, the application 
fip'-y 1 -* f ip(x)K(dx, y) is Lebesgue measurable. 
Mass conservation and physical interpretation lead to suppose 

n(dx,y) = l, JxK(dx,y) = ^, (21) 

so K,(y, .) is a probability measure and S Lf^ c (W + ). We moreover assume that the second moment 
of k is uniformly less than the first one 

r x 2 

/ -^K(x,y)dx <c< 1/2. (22) 

J y 

For the polymerization and fragmentation rates r and /3, we introduce the set 

V ■= {/ > : 3/x > 0, lim supx ^f{x) < oo and liminf x^f(x) > 0}. 

We consider 

/9 G n ^> 3r o > s.t. r E Lf oc (R+, x ro (i2;) n P (23) 

satisfying 

V-fT compact of (0, oo), 3m^- > s.i. t(x),(3(x) > for a.e. x £ K (24) 



3C>0, 7o>0 si. / y) dz < min 1,C - and — -SL^ 

Jo ^ ^V' ' T \ x ) 



(25) 



Notice that if Assumption (|25|) is satisfied for 70 > 0, then Assumption (I22p is automatically fulfilled 
(see Appendix A in |16j). We assume that growth dominates fragmentation close to L = in the 
following sense: 

~T ^ -^0 := {/> 3a > 0, / € L 1 (0,a)}. (26) 
We assume that fragmentation dominates growth close to L = +00 in the following sense: 

l im ^l = +00 . (27) 

Under these assumptions, we have existence and uniqueness of a solution to the first eigenvalue prob- 
lem; let us recall this result (see [27) [16]). 

Theorem [16J. Under Assumptions (|2ip -(j27 p . for < a, a < 00, there exist unique solutions (\,U), 
respectively to the eigenproblems ([2]), (J3J) and flD), and we have 

A > 0, 

x r rZY G L P (M+), Vr > - 7o , Vp G [l,oo], 
x r rW e W 1 ' 1 (R+), Vr > 0. 
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We also recall the following corollary (first proved in [27]). We shall use it at some step of our 
blow-up analysis. 

Corollary |27L 116] . Let r > and (3 > two given constants, 7, v G R such that 1 + 7 — > 0, 
and K£ verifying Assumptions (|21|) . (I22p and (|25 j) . T/ien t/iere exists a unique (#00)^00) solution to 
the eigenproblem ([2]) u>ii/i t(x) = tx u and f3(x) = fix 1 : 

Indeed, in this particular case, assumptions of the above existence theorem are immediate, and as 
already said both assumptions (|26|) and (|27p are verified if and only if 1 + 7 — i^>0. 



1.3 Proof of Theorem [2] 

It is straightforward to prove that K a satisfies Assumptions (j2ip -(l22l) and (|25h with the same constants 
c and C as k, thus independent of a. 
We have local uniform convergences in Ml r a — > tx v and (3 a — > fix 1 , if Assumption ([5]) holds, 

a— >L a—tL 

for L = as well as L = 00. 

For the two cases, the result is based on uniform estimates on v a and 9 a independent of a, in 
the same spirit as in the proof of the existence theorem (see [16]). When these estimates lead to 
compactness in L 1 (]R + ), we shall extract a converging subsequence, which will be a weak solution of 
Equation (I20p ; the global convergence result will then be a consequence of the uniqueness of a solution 
to Equation (|20p . This type of proofs is classical, and we refer for instance to [31] 116] for more details. 

As in [16], the proofs of the estimates that are needed to obtain compactness in L 1 are strongly 
based on Assumptions fl23]), ([25]), ([26]) and fl27]). Lemmas E] [3] and H prove respectively that j3 a , t q 
and K a satisfy Assumption (|27p . ([23]) and (|25p - ([26p uniformly for all a. 

Here we assume slight generalizations of Assumption ([5]), as precised in each lemma, in order to 
make appear more clearly where each part of Assumption ([5]) is necessary. 

Lemma 2. Suppose that 

3 f, 7 6 R s.t. l + 7 -^>0, and ^\ = 0(x^ 7 ) . (28) 

Then for all r > 0, there exist A r > and Ml a neighbourhood of L such that 

for a.e. x > A r and for all a E Ml, x @ a ( x ^ > r _ (29) 

T a (x) 

Proof. Let r > 0. 

For all a > 0, we define, thanks to Assumption (|27|). 



< p : — — - > r, for a.e. x > p > 

I ar(x) J 



p a := inf jp : 

Observe that p a is nondecreasing. Let e > 0, by definition of p a , there exists a sequence {£ a } with 
values in [1 — e, 1] such that 

77 x— < r - ( 30 ) 
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7" 

Thanks to Assumptions (1231) and (|24p . we have that — - G L?° (Ml) so that p a > +00 (else, 

xjj a^+oo 

since it is nondecreasing, it would tend to a finite limit, what is absurd by definition of p a ), and also 

T 

that — - > on Rl so that p n > 0. Hence, for some constant C : 

x/3 + a^O 

T J^\ <c£ aPa y^. (si) 

Then (|30p - (|3ip lead to, for some absolute constant C : 



'\ — £\ P a < ^ a P a <; 



a k a k 



1 k 



< C k r k 



which implies that lim sup —r < C k r k is finite. So we can define A r by 

A r := 1 + lim sup 

a— >L a 

Then for any x > A r we have, when a — > L, that a k x > and so 

a k xf3(a k x) ^ 
ar(a fc x) ~~ ' 

and by definition of f3 a and r Q it gives us the desired result. □ 
Lemma 3. Suppose that 

3v£R s.t. t(x) = 0(x u ) and 3r >0 s.t. x T H{x) G L,~(R+). (32) 

TTien /or a// A > and r > max (ro, — f), there exist C > and A/l a neighbourhood of L such that 

for a.e. x G [0,^4] and for all a G Ml, x r T a (x) < C. 

Suppose that 

3v£R s.t. r~ l {x) = 0{x~ v ) and 3 / u>0 s.t. inf x^t(x) > 0. (33) 

x^tL x£[l,+oo) 

for all e > and m > max(/x, — 1/), i/iere exisi c > and Ml a neighbourhood of L such that 

for a.e. x > e and for all a G Ml-, x m T a {x) > c. 

Proof. We have to treat separately the case L = and L = +00. 

Let us start with L = 0. Notice that in this case, if t(x) = 0(x u ) then thanks to Assumption (|23p 
ro = — 1/ suits to have (132p . Considering r > -y, we have for some constant C > : 



sup ess (x r Tc,(x)) = sup ess (a k ^ r+v ^y r T(y)) 

xG[0,A] y£[0,a k A] 

< C sup (a- k{r+u) y r+u ) = CA r+u . 

a~>0 [0,a k A] 
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For m > max(/i, — v) and using Assumption (j33|) we have for some constants ci,c 2 > : 
inf ess (x m r a (x)) = infess {a- k[m+u) y m T{y)) 

x£[e,oo) ye[a fe e,oo) 

> min (infess (a- k{m+u) y m T(y)), infess {a~ k{:m+v) y^T{y)) 

\[a k s,l] [l,oo) 



fci inf (a- fc ( m+ ^y m+ ^),c 2 Q- fc(m+i/) 
V [o fe e,l] 



> min(ci inf {a~^ m ^>y m ^),c 2 a-^' ,L ^> )= mm ( Cl e m ^ ,c 2 a 

Now we consider L = +00 and r > max(ro,— v). Thanks to Assumption (|32|) we have for some 
constants C\ , C2 > : 

supess (x r r Q (x)) = sup ess (a~ k ^ r+u ^y r T(y)) 

x£[0,A] y£[0,a k A] 

< sup ess {a- k[r+u) y ro T{y)) + sup ess (a- fc(r+l Vr(y)) 

[0,1] [l.ctM] 

< Cia- fc(r+l/) + C 2 A r+u . 

For m> —v and using Assumption (I33D . we have for some c > : 

infess (x m T a {x)) = infess (a" fc(m+1 VV(n)) 

x£[e,oo) ye[a k e,oo) 

> c inf { a - k{m+y) y m+u ) = ce m+u . 

a— >-oo [a k e,oo) 



Lemma 4. Suppose that 



lo r a (x) 

Proof. Thanks to Assumption [M] we have for some constant C > : 



fc - - "£Q fe 



dx = - '^-+dx < C- x i ~ v dx = Cke k . 



1 



□ 



0(x) 

3 u, 7 e M s.fc 7 + 1 - > 0, and -f^ = 0(^ 7_v )- (34) 
TTien /or a// p > 0, i/iere exisi e > and A/l a neighbourhood of L such that 

Va G Ml, / ^r^<p. 



o Tq,(x) a J t(x) a-^L a J 

The result follows for e small enough. □ 



Thanks to these preliminary lemmas, we can prove Theorem [2j 

Proof of Theorem^ We prove estimates on {v a } and {6 a }, uniform in a — > L, in order to have 
compactness. Then the convergence of the sequences is a consequence of the uniqueness of (600, v oo)- 
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tion (122 



First estimate: Abound for x r v a , r > 0. For r > 2, we have by definition and thanks to Assump- 

ry x r ry x 2 

/ —K a (x,y)dx < / -^K a (x,y)dx 
Jo V Jo U 



f-a k y x 2 



= / 7 k , 2 K{x,a k y)dx < c. 

Jo {a k yr 

So, multiplying the equation (|1T|) on v a by x r and then integrating on [0, oo), we find 
/ (1 — 2c)x r f3 a (x)v a (x) dx < r x r ~~ 1 T a (x)v a (x) dx 



= r I x r l T a {x)v a (x) dx + r I x r T a (x)v a (x) dx. (35) 

Jx<A Jx>A 

Thanks to Lemma [2j and choosing A = Az. in (|35p with u < 1 — 2c, we obtain 

'sup^,)!^- 1 ^} 



■ 

x r (3 a (x)v a (x) dx < 



1 - 2c - u 

and the right hand side is uniformly bounded for r — 1 > max(ro, —v) when a — > L, thanks to LemmaH 
Finally 

Vr > max(2, 1 + r , 1 — u), 3C r , J x r f3 a (x)v a (x) dx < C r . (36) 

Moreover for all a we have f v a dx = 1. So, using Lemma with /3~ 1 (x) = 0{x~^) instead of 
t~ 1 {x) = 0(x~ u ), we conclude that uniformly in a — > L 

x r v a £ L 1 (R+), Vr > 0. (37) 



Second estimate: 9 a upper bound. The next step is to prove the same estimate as (|36p for < r < 
max(2, 1 + ro, 1 — v) and for this we first give a bound on T a v a . Let m = max (2, 1 + ro, 1 — v), then, 
using e and p < | defined in Lemma d] and integrating (|1T|) between and x < e, we find (noticing 
that the quantity T a (x)v a (x) is well defined because Theorem [16] ensures that r a v a is continuous) 



T a (x)v a (x) < 2 J J P a (y)v a (y)K a (z,y) dy dz 

< 2 J (3 a (y)v a (y) dy 

pa rco 

= 2 J P a (y)v a (y)dy + 2 y Pa{y)v a (y)dy 

< 2snp{r a v a } f ^ dy + 2e~ m T 

(0,e) JO T a{y) Jo 



(0,e) 

< 2p sup{r Q f Q } + 2e _m C„ 

(0,e) 



y m p a (y)v a (y)dy 



Consequently, we obtain 



1 + 2C e~ m 

sup r Q! (a;)'y C( (x) < — ^ := C. (38) 

xG(0,e) 1 _ l P 
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Then we can write for any < r < m 

x r f3 a (x)v a (x) dx = / x r f3 a (x)v a (x) dx + / x r (3 a (x)v a (x) dx 



< e r snp{r a v a } dx + e r ~ m [°° x m ^ a (x)v a {x) dx 

<0,e) Jo T a[ x ) J e 



(0,e) 

< Cpe r + C m e r - m :=C r . 

Finally we have that 

Vr>0, 3C r , J x r (3 a (x)v a (x) dx < C r (39) 

and so 

6 a = j f3 a v a < C . (40) 

Third estimate: L°° bound for x~ l0 T a v a . First, integrating equation (|17p between and x we find 
r a (x)v a (x) < 2 / p a (y)v a (y) dy = 26 a < 2C , Vx > 0. (41) 



It remains to prove that bounded in a neighborhood of zero. 

Let us define f a : x h- > sup( 0jX ) T a v a . If we integrate ([2|) between and x' < x, we find 

T a (x')v a (x') <2 J J j3 a (y)v a (y)K a (z,y)dydz<2 j j f3 a (y)v a (y)K a (z,y) dy dz 
and so for all x 

f a (x) < 2 / / /3 a (y)v a (y)K a (z,y) dy dz. 



'0 



Considering e and p from Lemma 0] and using (|25p , we have for all x < e 
f a (x) < 2 j j f3 a (y)v a (y)K a (z,y) dy dz 
= 2 [ (3 a (y)v a (y) J K a (z,y)dzdy 

P a (y)v a (y)min(l,C^-y ) dy 
2/ P a {y)v a (y)dy + 2C P a (y)v a (y)(-) dy + 2C P a (.y)v Q (y)(-) dy 

Jo Jx K V J Je K V J 

2 



< 2 



T a {y)v a (y)dy + 2Cx^° — — dy + 2C P a {y)v a {y)[-) dy 

J x T a {y) y J0 Je K y J 



lo My) 

< 2f a (x) [ e !!M d y + 2Cz» f ^MhM d y + 2Ce^WaV ct \\ L ,x 

Jo MV) Jx T a{y) V 

If we set V a (x) = x~ 10 f a (x), we obtain when a — > L 

(1 - 2p)V a (x) <K £ + 2C f ^44v Q (y) dy 

Jx MV) 
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2Cp 

K e 1 ~ 2p 

and, thanks to the GrdnwalFs lemma, we find that V a {x) < . Finally 

1 — 2p 

2Cp 

sup{x-^T a (x)v a (x)} < Keel ~ 2P . (42) 

(0,e) 1 - 2 P 

The bound (|41|) with Assumption (|24|) and the bound (|42p with Lemma[H(in which we replace /3 a ( x ) by 
x 70 ) ensure that the family {v a } is uniformly integrable. This result, associated to the first estimate, 
gives that {v a } belongs to a compact set in L 1 -weak thanks to the Dunford-Pettis theorem. The 
sequence {0 a } belongs also to a compact interval of R+, so there is a subsequence of {(v a , 6 a )} which 
converges m L x -weak x R. The limit is a solution to (I20j) , but such a solution is unique so the sequence 
converges entirely. To have convergence in L 1 -strong, we need one more estimate. 

Fourth estimate: W 1,1 bound for x r r a v a , r > 0. First, the estimates (|37p and (|42p ensure that x r r a v a 
is uniformly bounded in L 1 for any r > —I. Then, equation (|17|) gives immediately that 

— (x r T a (x)v a (x)) | dx < r J x r ~ 1 T a (x)v a (x) dx + 9 a + 3 J I3 a (x)v a (x) (43) 

is also uniformly bounded. For r = 0, the same computation works and finally bounded in 

W 1 ' 1 (R + ) for any r > 0. 

The consequence is, thanks to the Rellich-Kondrachov theorem, that {x r T a v a } is compact in L 1 - 
strong and so converges strongly to tx t ' +u v 00 (x). Then, using Lemma [3] and estimate (j4"2j) . we can 
write 



x r \v a (x) — Voo(x)\ dx < / x r \v a (x) — Voo(x)\ dx + / x r \v a (x) — Voo(x)\ dx 







6 x r+7 ° Z" 00 



< C I — -— dx + C I x r \T a (x)v a (x) — tx^VoqIxM dx. 

J0 T a (x) 



E 



The first term is small for e small thanks to Lemma [4] and the second term is small for a close to L 
due to the strong L 1 convergence of {x r T a (x)v a (x)} . This proves the strong convergence of {x r v a (x)} 
and ends the proof of Theorem [2j 

□ 



2 Further Results 

2.1 Critical case 

In the case when lirn r _ > o (3(x) or lim x _5.o is a positive constant, we can enhance the result of 
Theorem [T] if we know the higher order term in the series expansion. Assumptions (|44|) and (|46p of 
Corollary Q] are stronger than Assumption ([5]) , but provide a more precise result on the asymptotic 
behavior of X a , A„. 
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Corollary 1. If (3 admits an expansion of the form 



P(x) = /3 + Pi* 11 + o (x 71 ), 71 > (44) 

x— >0 



with /3q > and f3\ ^ 0, then we have the following expansion for A 



K = A) + ( Pi I x^ Voo {x) dx ) a k ^ + o (a fe71 ). (45) 

a— >0 



In the same way, if r admits an expansion of the form 



rlx) = TqX + T X X Vi + o (x Ul ), Vi > 1 (46) 
x-^-0 



with tq > and T\ ^ 0, then we have 



K = ro + ( Tl lf^±) ^-D + (a ^-D). (47 ) 



m-oo 



Proof. First we assume that /3 admits an expansion of the form (|44p and we want to prove (|45|) . By 
integrating Equation (|17p we know that A Q f = f f3 a (x)v a (x) dx and so multiplying by a _fc71 
gives 

a-*^(Aa - A>) = / «- fe71 (/3a(x) " A))«a(x) dx. 



Now the proof is complete if we prove the convergence 

a~ k ^((3 a (x) - po)v a (x) dx > [ pix^v^ix) dx. (48) 

Q^O J 

For this we use the expansion (j44j) which provides, for all x > 0, the following one 

/S«(aO = /?o + /9ix 7l a fc71 +o(a fc71 ). (49) 

a— i-0 

Let m > 7i such that limsup:r _m /3(x) < oo (see Assumption ()23p ) and define 

x— >oo 

a- k ^(p a (x)-/3 ) 



f(% '. x i y 



3.71 _|_ x m 



We know thanks to (|49p that f a (%) — > for all x. Moreover we have thanks to Theorem [2] 

a^O xfi -f x m 

that (x 71 + x m )v a (x) > (x 71 + x m )v 00 (x) in L . So we simply have to prove that f a is uniformly 

a— i-0 

bounded to get (I48p (see Section 5.2 in [21] )■ Thanks to (I44p and the fact that limsupx _m /3(x) < oo 

X— >OD 

with m > 7i > 0, we know that there exists a constant C such that 

|/3(y)-/3 |<C(y 71 +y m ), Vy > 0, 

and so, because a — > 

a- fe7l |/3(y) - A)| < C(a- fc7l y 71 + a~ km y m ), 

what gives, for x = a k y 

a- k ^\P(a k x) - A) | < C(x 71 +x m ) 
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and we have proved that f a (x) < C. 

The same method allows to prove the result on A , starting from the identity 

(A a - r ) J xw a (x) dx = J (T a (x) - T x)w a (x) dx 

and using the fact that (146ft provides the expansion 

r a (x) = r x + Tix i V^ 1 - 1) +o(a' (l/1 - 1} ). 

a— >oo 

□ 

2.2 Relaxed case 

Here we relax Assumption ([5]) and look if the asymptotic behavior of X a and A a obtained in Theorem [1] 
remain true. The case we are the most interested in is the case when the limits are zero (see the 

t(x) 

applications at Section[3|): is the condition lim (3{x) = (resp. lim = 0) necessary and sufficient 

x^tL x 

to have lim A Q = (resp. lim A n = 0)? The following theorem gives partial results in the direction 

a— ¥L a— >L 

of a positive answer to this question. The assumptions required are weaker, but the results also are 
less strong: we obtain asymptotic behavior for the eigenvalue but cannot say anything yet on the 
eigenvector behavior. 

Theorem 3. Let us suppose that all assumptions of Theorem^ are verified except Assumption ([5]). 

1. If t(x) = o(x u ), P = <9(x 7 ) and f3(x)~ 1 = 0(x~ 7 ) with 7 + 1 — v > 0, we have that 

x— s-0 x— >0 x— >0 

7 

if'Y > 0, then lim X a = 0, and more precisely X a = o(a 1 +~t-"), 

a— s-0 «— s-0 

l-u 

if v > 1, then lim A n = 0, and more precisely A a = o(a 1+ i- v ). 

a— >-oo a— s-oo 

2. If (3(x) = o(x~f) and r(x)" 1 = 0(x~ u ) with 7 + 1 - v > and 7 < (so that v < I), we 

x— yoo x— >oo 

have that 

7 \ — v 

lim Aq, = lim A a = 0, and more precisely \ a = o(a 1+ ~<- v ), A a = o(a 1+ "'- v ). 

a— s-oo a— >0 a— >+oo a— >0 



Remark. In the second assertion of this theorem, we notice that Assumption (j27|) means = 

X X— YOO 



t{x) 

. ' ' . so the condition f3(x) = o(x 7 ) with 7 < imposes lirn^oo = 0. 
Proof of Theorem{3$Jl We perform the dilation defined by ©: v a (x) = a k U a (a k x) with k - 



Thanks to the assumption /3(x) _1 = 0(x -7 ) and r = 0(x u ), the conclusions of Lemma [2] and 

z-s-0 

Lemma [3] still hold true. Hence we have the following bound (see the first estimate in the proof of 
Theorem [2]) 

"up^^af-Va) 
x r p a (x)v a {x) dx < - ^ - 
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where T a (x) and (5 a are defined by (|13p . and the right hand side is bounded uniformly in a for 
r > max (2, 1 + tq,1 — u). Let e > and write 



a i+7-" A a = a = / /3 a Va < / &*(a;)?; a (a;) ds + e r y x r /3 a (x)v a (x) dx 

/sup( 0)Aj . ) (x r-1 Ta) 



< sup f3 a + e 

Thus, since sup £ ro T a — >• 0, we obtain 



(o,e) l-2c-u 



(0,4r ) 



a->0 



hmsup# Q < lim sup sup /3 Q < Ce 7 . 

a->0 (0,e) 

This is true for all e > so the assertion Q] of Theorem [3] is proved (the same proof works with the 
fragmentation parameter a). 

□ 

Proof of Theorem\MM We perform the dilation defined by (jBJ) v a (x) = a k U a (a k x) with k = 
Thanks to Assumption ()34|) for L = +oo, we still have the conclusion of Lemma [Hand it is sufficient 
to bound r a v a on (0,e) for e > 0. We refer to the proof of Theorem [2] in Section [L3| second estimate, 
and write: 



r a (x)v a (x) < 2sup{r a v a } [ ^T^r dy + 2 / 
(0,e) i0 



dy + 2 Pa(y)v a (y) dy 

(0,6) 

< 2p sup{r a v a } + 2 sup /3 a . 

(0,e) (e,+oo) 

Taking for instance e small enough so that p < \ in this estimate, and a large enough so that 
sup( £j+0O ) f3 a < C, we obtain the boundedness of T a v a on (0, eo). Then we write for e > 0, 



a 1 +7— A„ = (9 Q = / 



< sup r a v a I — + sup /3 a . 



(0,e) JO T a ( e ,oo) 



<c >o 

fc u a— ¥oo 

The latter estimate is a consequence of the assumption /3(x) = o(x 7 ). 
We do the same computations for the parameter a. 

□ 



2.3 Generalized case 

In this section, we completely free ourselves of Assumption ([5]) and give some results about the 
asymptotic behavior of the first eigenvalues for general coefficients. The technics used are completely 
different than the self-similar ones. But the results still give comparisons between A Q and /3(x), and 
between A a and 
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Theorem 4. 1. Polymerization dependency. 

7~(x) 

If f3 G L^ C (M+) and limsup < oo, then lim sup^^g X a < limsup^o f3(x). (50) 

X— ¥OD X 

If - G L\ :={/, 3a > 0, / G L^O,^}, i/ien liminf^ooAa > liminf^oo /3(x). (51) 
2. Fragmentation dependency. 

Tl X I 

If (3 £ L^ C (]R + ), i/ien there exists r > suc/i i/iai lim sup A a < r lim sup . (52) 

a— >-0 x— >oo 3? 
// liminf /3(x) > 0, and if — G Lq, i/ien lim A a = +oo. (53) 

x— >oo T cn-oo 

We first state a lemma which links the moments of the eigenvector, the eigenvalue and the polymer- 
ization rate. 

Lemma 5. Let (U,X) solution to the eigenproblem (|2|). For any r > we have 

J x r U{x) dx<jj x r ~ 1 T(x)U(x) dx. (54) 
Proof. Integrating Equation ([2]) against x r we find 

rx T ~ l T(x)U(x) dx + A / x r U(x)dx+ \ x r j3(x)U{x) dx 



x 



2 / x r / p(y)n{x,y)U(y)dydx 



n 



f) 



2 / /%)W(y) / x r K(x,y)dxdy 



o 



< 2 J P{y)U(y)y r 1 / xn(x,y)dxdy 
y r f3(y)U(y)dy. 



o 



□ 

Proof of Theorem [Tpl ([50] ) . We only have to consider the case limsup^Q P( x ) < °°- 111 this case, 
/? G L^ C (K+) since it is assumed that f3 G L^ C (M+). So for e > 0, we can define j3 e := sup( 0)£ ) f3(x) < oo 
and, thanks to Assumption (|23p . there exist positive constants C £ and r such that /3(a;) < Cx r for 
almost every x > e. As a consequence, we have by integration of Equation ([3]) 



\ a = (3{x)U a (x) dx < P £ + C / x r U a (x) dx. 



We can consider that r > ro + 1, with ro defined in Assumption (j23D . and thus Lemma [5] and 
Assumption lim sup,^^ < oo give 

y x r U a {x)dx <Y J xr ~ lr ( x Pa(x) dx < y C U + y £74(x)dx^ 
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for a new constant C. Combining these two inequalities we obtain 

K < aC ( 1 + ; , 1 <aC[l 



f x r U a (x) dx J \ A a — /3< 

Then, either \ a < /3 E , or we obtain by multiplication by A Q — j3 £ > that 

\ 2 a - (We + «C)A a - (1 - ^)aC < 0, 

and so 

A„, < 

2 

Finally we have 

lim sup A a < /3 £ 

and this is true for any e > 0, so 

lim sup A a < limsup/3(x). 

a->0 x-»0 



A Q < - ( & + aC + V /3 e 2 + 4aC + a 2 C 2 



□ 



Proof of Theorem \4$H (EU) • Let A > and define /3a := hif (a,oo) P- Since - € Lq and thanks to 



Assumption (|24p we can define /a := Jo* hIy ^ 00 • Then we have by integration of Equation (|3|) 



Aa = J (3(y)U a (y)dy 

poo 

> 13a / W a (y)dy 



A 

A 



§a(i~J U a {y)dy). 



We know, by integration of Equation ([3]) between and x, that for all x > 0, a.T(x)U a (x) < 2X a . 
Thus we obtain 



A Q > A* 1- / 2A U 

ar{y) 



o 



which gives 

Aa ^ — o~~ 

" 1 + a-^L 7 

and 

liminf A a > /3a • 

a— >oo 

Finally, since the previous inequality is true for all A > 0, we find 

liminf Aq > liminf /3(x). 

a— >oo u— >+oo 

□ 
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Proof of Theorem ([52]) . The fact that /3 G Lf£ c (R + ) associated to Assumption ([23]) ensures the 
existence of two positive constants C and r such that for almost every x > 0, < C(l + a? r )- So, 

integrating Equation (jlj), we have 



K = aJ P{x)V a (x) dx < aC (l + J x r V a 



(x) dx I . 



To prove ([52]) . we only have to consider the case lim sup^^ < oo. So, for any A > 0, we can 

define Ta '■= su Px>A ~^r~ < 00 an< ^ we have, thanks to Lemma [5] and considering r > ro + 1 where ro 
is defined in Assumption (|23p . 

z r V„(x)da; < x r-1 r(s)V B (s) rfs < (^C + rJ ^ x r V a (x)dx 

Combining the two inequalities we obtain 

( C \ ( aC 

A < r [ta + r / n , <r[TX + 



J x r V a {x) dx) ~ \ A a - aC 
Then, either A a < aC, or we obtain by multiplication by A a — aC > that 

A„ - (rfX + aC)A a - r(l - TX)aC < 0, 

and so 



Aa < \ (tta + oC + V(^) 2 + 4raC + a 2 C 2 ) 



In the two cases we obtain that 



lim sup A < rTA, 
a->0 



□ 



and, doing A — > oo, we conclude 

f( x) 

lim sup A a < r lim sup . 

a— >0 x— >ao X 



Proof of Theorem^ (153j) . Let e > 0. Since lim inf j.-^ f3(x) > 0, we have thanks to Assumption (1241) 
that j3 £ := inf( £)00 ) /5 > 0. Since ~ G Lq, we have thanks to Assumption (|2"I]) that I e := J* e < oo 
and that lim £ _>.o I £ = 0. 
By integration of equation @, we find 

A = a J P(y)V a (y)dy 

V a (y) dy 

= a§e(l - f v*{y)dy). 

We know, as previously by integration between and x, that for all x > 0, T(x)V a (x) < 2A a . Thus we 
obtain 

dy 



K>aP £ (l- / 2A a 
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which gives 

An > 



a " 1 + 2a&/ £ 
Finally we have 

lim inf A a > > +00 

a->oo 2I £ e->0 

and it ends the proof of Theorem SJ ([53)) . □ 



3 Applications 

Before giving applications, we remind a regularity result whose proof can be found in |28j . 

Lemma 6. Under the assumptions of Section \1 . H the functions Qh->A a and a 1— > A a are well defined 
and differentiate on (0, 00). 

3.1 Numerical scheme based on Theorem [2] 

Let us first present the method we use to compute numerically the principal eigenvector A without 
considering any dependency on parameters. Then we explain how the self-similar change of variable © 
and the convergence result of Theorem[2]can be used to compute the dependencies a 1— > A Q and a 1— > A a , 
when parameters a and are very large or very small. 



The method used to compute A solution to Equation ([2]) is first to compute a numerical approxi- 
mation of the first eigenvector U, and then use the identity 



oo 



A= / P(x)U(x)dx. 









The General Relative Entropy (GRE) introduced by [29j [30j EI] provides the long time asymptotic 
behavior of any solution to the fragmentation-drift equation ([I]). At sufficiently large times, these 
solutions look like U(x)e xt where U and A are the eigenelements defined at ([2]). More precisely we 
have 

fOO 

\u(x,t)e~ xt - {u(-,t = 0),4>)U(x)\4>(x)dx > 0, 

t— >oo 

where <p is the dual eigenvector of Equation ([2]) (see [TUEO] for more details) and (u, 0) = u(x)4>(x) dx. 
It is even proved in \24\ IS], under some assumptions on the coefficients, that this convergence occurs 
exponentially fast. We use this convergence to compute numerically the eigenvector IA. We consider, for 
uo £ L 1 (]R + ) an initial function satisfying f^°uo(x)dx = 1, the solution u(x,t) to the fragmentation- 
drift equation (TTJ. Since we do not know yet the value of A, we define the normalized function 

-/ n, u(x,t) 
u{x,t) := -^5o - — . 

J u(x,t)dx 

We can easily check that u satisfies the equation 







d t u{x,t) + d x (r(x)u(x,t)) + ( / (3(y)u(y,t) dy ) u(x,t) + /3(x)u(x,t) = 2 \ f3(y)K(x,y)u(y,t) dy, 

(55) 
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with the boundary condition r(0)u(0, t) = 0, and that the convergence occurs 



\u(x, t) - U{x)\(j){x) dx > 0. (56) 



The scheme used to compute IA is based on the resolution of Equation (|55p for large times and the use 
of (|56p for the stop condition. 

Numerically, Equation (|55|) is solved on a truncated domain [0, R] so the integration bounds have 
to be changed and we obtain, for x £ [0, R], 

d t u(x,t) + d x (r(x)u(x,t)) + (J P(y)u(y,t) dy\ u(x,t) + P{x)u(x,t) = 2 J (3(y)n(x,y)u(y,t) dy. 

(57) 

What we loose when we solve this truncated equation are the integral terms f R j3(x)u(x, t) dx and 
J R (3(y)n(x,y)u(y,t) dy, and the outgoing flux T(R)u(R,t) at the boundary x = R. To be as close as 
possible to the non-truncated solution, we have to choose R large enough so that these quantities are 
small enough. It is proved in [16] that f3{x)U{x) and t{x)U{x) are fast decreasing when x — > +oo. 
The value of R has to be adapted to have, when u is close to the equilibrium U, values of t(x)u(x, t) 
and (3(x)u(x,t) smaller than a fixed parameter e for x close to R. Parameter e is expected to be very 
small, it is also used for the stop condition (|59p . 

We assume that [0, R] is divided into N uniform cells and we denote Xi = iAx for < i < N with 
Arc = ^. The time is discretized with the time step At and we denote t n = nAt for n € N. We adopt 
the finite difference point of view, namely we compute an approximation uf of u(xi, t n ). It remains to 
explain how we go from the time t n to the time t n+l . To enforce that = 1 at each time step, 

we split the evolution into two steps. First we compute, from (uf )i<i<N, a vector (ti™ +1 )i<j<7v which 
is obtain with the formula 

..n+1 ~n -n+1 _ ~n N 

= ~ /T 1 + 2Ax £ 3 jKi:j ^ (58) 

where /3j = /3(xj), Tj = r(xi) and m ~ = K>(xi,Xj). This is a semi-implicit Euler discretization of the 
growth- fragmentation equation ([1]). We choose this scheme to ensure the stability without any CFL 
condition, since the scheme is positive. Then we set 

,.n+l 



' Ax EyLi u ] +1 

and we have that the discrete integral of (u" +1 )i<i<Ar satisfies Ax YaLi = Ax J2iLi^ L i' = 1- 
Inspired from the L 1 convergence (|56|) . we end the algorithm when 

%Y.K~ur l \<e (59) 
i=i 



where e <C Ax. Then we have 



N 



A~ AxJ2^i- 



i=l 
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The semi-implicit scheme (|58p is efficient to avoid oscillations on the numerical solution but it is not 
conservative. It has to be avoided if we want to solve Equation ([I]) for any time. Here we are only 
interested in the steady state of Equation (|55|) , so the non conservativity does not matter because the 
steady state is the same for an implicit or an explicit scheme. 

Now we want to compute \ a and A for a large range of a and a. According to the discussion 
in introduction, the eigenvectors U a and V a are concentrated at the origin for a small or o large 
and, conversely, it is spread for a large or a small. Then, to avoid an adaptation of the truncature 
parameter R or an adaptation of the discretization size step Ax when ex and a vary, we compute a 
defined in Equation (|17p . To this end, we have to compute the dilated eigenvector v a defined in ([6]) 
which converges to a fixed profile Voo when a — > L, as stated in Theorem [2j This convergence ensures 
that the vector v a does neither disperse nor concentrate too much when a varies, and so we can find 
a truncature and a size step which suit for any a — > L. It remains to distinguish L = from L = +00 
by dividing (0, +00) into two sets: for instance (0, 1] and (1, +00). For < a < 1, we use the dilation 
coefficient k associated to v and 7 such that r(x) ~ x v and /3(x) ~ x 7 . For a > 1 we do the 

x— >0 x— 5-0 

dilation associated to v and 7 such that t(x) ~ x v and /3(x) ~ x 7 . Finally, we use Equation (|19p 

x— >oo x— >oo 

in Lemma Q] to recover \ ct or A a from the numerical value of a . 



All the figures of the paper are obtained with this numerical scheme. 



3.2 Steady States of the Prion Equation 

To model polymerization processes, Equation (pQ) can be coupled to an ODE which leads the evolution 
of the quantity of monomers. It is the case in the so-called "prion equation" (see [22| [3i]) 



poo 

/ r(x)u(x,t) 
Jo 



dx 



J^(x, t) = -V{t)-^-(T(x)u(x, £)) - \p(x) + n(x)]u(x, t) + 2 J x £(l/)«0c, y) u(y, t) dy, 
u(0,t) =0, 

where the quantity of monomers is denoted by V(t). In this model, the monomers are prion proteins, 
produced and degraded by the cells with rates £ and 5, and attached to polymers of size x with 
respect to the rate r(x). The polymers are fibrils of misfolded pathogenic proteins, which have the 
ability to transconform normal proteins (monomers) into abnormal ones by a polymerization process 
not yet very well understood. The size distribution of polymers u(x, t) is solution to the growth- 
fragmentation equation ([T]) in which V(t) is added as a multiplier for the polymerization rate. A 
degradation rate (jl{x) is also considered for the polymers. For the sake of simplicity, this rate is 
assumed to be size-independent in the following study (/u(x) = fj,o). 

Equation (|60p models the proliferation of prion disease. An individual is said to be infected by prion 
disease when polymers of misfolded proteins are present, namely when u(-,t) ^ at the present time 
t. 

The coupling between V(t) and u(t, x) appears in the equation for u as a modulation of the poly- 
merization rate. One sees immediately the link with the eigenproblem ([3]) satisfied by U a . Indeed, U a 
is the principal eigenvector linked to the linearization of the prion equation around a given monomer 
quantity V = a. Investigating the dependency of the fitness Ay with respect to the polymerization and 
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fragmentation coefficients is a first step towards a better understanding of the disease propagation. 
Indeed it has been reported that the course of prion infection in the brain follows heterogeneous pat- 
terns. It has been postulated that the neuropathology of prion infection could be related to different 
kinetics in different compartments of the brain [12] . 

Modelling the propagation of prion in the brain requires a good understanding of possible dynamics 
(e.g. monostable, bistable etc). Such a study can be done through the dependency of the first 
eigenvalue on parameters. In [8j , it is shown that, under some conditions, the coexistence of two 
stable steady states can happen (one endemic and one disease-free). A steady state (VocU^x)) is a 
solution to 







S + I t(x)u oc (x) dx 
o 







fi Uoo(x) = -Voo — ^^u^x)) - P(x)uoo(x) + 2 / /3(y)K(x,y)u 00 (y)dy, 



(61) 



ltoo(0) 



0. 



The disease-free steady state corresponds to the solution without any polymer \ V = |,u = 0! 

Other steady states can exist and are called endemic or disease steady states. They are solution 
to System (jSTl) with Voo > and Uoo ^ nonnegative. To know if such disease steady states exist, 
we recall briefly here the method of [HJ [9]. A positive steady state Mqo can be seen as an eigenvector 
solution of (El) with a = VJx, such that 



A, 



Ay. 



(62) 



It shows the crucial importance of a study of the map V i- > Xy . Any value Voo solution to ([62]) provides 
a size distribution of polymers 

Uoo(x) = QooU Voo (x). 

The quantity of polymers Qoo is then prescribed by the equation on monomers and has to satisfy the 
relation 

e^ 1 - * 

°° Voo J t(x)Uv 00 (x) dx' 
The quantity of polymers has to be positive, what is equivalent to the condition 



Vr 



« < | = v. 

o 



Finally, the disease steady states correspond exactly to the zeros of the map V i- > Ay — fio in the 
interval (0, V). Thanks to the different results of Sections Q] and [21 we know that this map is not 
necessarily monotonic, as assumed in [9]. Thus, by continuity of the dependency of A on V (see 
Lemma [6]), there can exist several disease steady states for hq well chosen and V = I large enough. 
This point is illustrated on the example below, for which there exist two disease steady states. 

We can investigate the stability of the disease-free steady state through the results obtained in [HI [9]. 
For this, we introduce the dual eigenvector Tp of the growth-fragmentation operator with transport 
term V 

— d f x 

-Vt(x)—(Tp(x)) + (J3(x) + X)7p(x) = 2(3(x) / n(y, x)7p(y)dy, x > 0, 

Jo (63) 



<p(x) > 0, J °° (p(x)U v (x)dx 



1. 
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We assume that we are in a case when there exist two constants K\ and K2 such that 



< K±(p(x), and t(x) < K2 i p{x). 



(64) 



This assumption generally holds true when is bounded because (p grows linearly at infinity ac- 
cording to general properties proved in [TU ETJ EU [32] . Then we can reformulate the theorems of JHl [9] . 

Theorem [9] (Local stability). Suppose that assumption (|64p holds true and that Xy < fiQ. Then the 
steady state (V, 0) is locally non-linearly stable. 

Theorem [8] (Persistence). Suppose that assumption (f64"|) holds true, V(0) < V, / °°(1 + x)u(t, x) dx 
is uniformly bounded, and that Xy > (J,q. Then the system remains away from the steady state (1^,0). 
More precisely we have: 



lim inf 

t— >oo 



ip{x)u{x, t) dx > 0. 



Example. Let us consider the same coefficients as in Figure 2(a) We can choose /Uq small enough 



to ensure the existence of two values V\ < V<i such that Ay = Ay, = As a consequence, we know 
thanks to the previous study that there exists no disease steady state if V < V±, one if V\ < V < V2, 
and two if V > Vi- Concerning the stability of the disease- free steady state, we first notice that 
the fragmentation rate /3(x) satisfies the assumption f3(x) < A + Bx which is sufficient to have that 
J °°(l + x)u(t,x) dx is uniformly bounded (see [8] Theorem 2.1). Thus we can apply the previous 
theorems and we have that (V, 0) is stable if V < V±, unstable if V± < V < V2, and recovers its (local) 
stability if V > V2. In Figure EJ the graph of the negative fitness V 1— > fj,o — Ay is plotted (because the 
quantity of polymers influences the evolution of V{t) with a negative contribution) and the zones of 
stability and unstability for V are pointed out. The non intuitive conclusion is that an increase of the 
production rate £ or a decrease of the death rate 5 can stabilize the disease-free steady state. What 
happens in this situation is that the largest polymers are the most stable since hm a ;_5. 00 /3(x) = (this 
situation is biologically relevant, see for instance [35]). When the quantity of polymers is large, the 
polymerization is strong and it forms such long stable polymers. Because they do not break easily, 
their number does not increase very fast, i.e. the fitness of the polymerization-fragmentation equation 
is small. But the degradation term is assumed to be size-independent, and then the fitness Ay becomes 
smaller than /xq for V large enough. This phenomenon stabilizes the disease-free steady state because, 
when polymers are injected in a cell, they tend to disappear immediately since A^ < 

Concerning the stability of the disease steady states, the study is much more complicated. Never- 
theless, we can imagine that V\ is stable and V2 unstable. This postulate is based on Figure [3] and on 
the results obtained in [22j [23] in a case when System ([60]) can be reduced to a system of ODEs. 



3.3 Optimization of the PMCA protocole 

Prion diseases, briefly described in Section [3.2l (see [25| for more details), are fatal, infectious and neu- 
rodegenerative diseases with a long incubation time. They include bovine spongiform encephalopathies 
in cattle, scrapie in sheep and Creutzfeldt- Jakob disease in human. It is then of importance to be able 
to diagnose infected individuals to avoid propagation of the disease in a population. But the dynamics 
of proliferation is slow and the amount of prion proteins is low at the begining of the disease. More- 
over, these proteins are concentrated in vital organs like brain, and are present in minute quantities in 
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Figure 3: The negative fitness V \— > fj,o — Ay is plotted for the same coefficients as in Figure 2(a) The 
zeros V\ and V2 correspond to disease steady states and separate the areas of stability or unstability 
of the disease-free steady state. 



tissues like blood. To be able to detect prions in these tissues, a solution is to amplify their quantity. 
A promising recent technique of amplification is the PMC A (Protein Misfolded Cyclic Amplification). 
Nevertheless, this protocole is not able to amplify prions for all the prion diseases from tissues with 
low infectivity. It remains to be improved and mathematical modelling and analysis can help to do 
so. 



PMCA is an in vitro cyclic process that quickly amplifies very small quantities of prion proteins 
present in a sample. In this sample, the pathogenic proteins (polymers) are put in presence of a large 
quantity of normal proteins (monomers). Then the protocole consists in the alternance of two phases: 

- a phase of incubation during which the sample is left to rest and the polymers can attach the 
monomers (increase of the size of polymers), 

- a phase of sonication during which waves are sent on the sample in order to break the polymers 
into numerous smaller ones (increase of the number of polymers). 

To model this process, we can use the growth- fragmentation equation ([I]) as in (|60p . The main differ- 
ence is that the PMCA takes place in vitro, then there is no production of monomers. The monomers 
are in large excess to improve the polymerization, so we can assume that their concentration is constant 
during the PMCA. It remains to introduce the "sonication" in the equation. Because the sonication 
phase increases the fragmentation of polymers, a first modelling can be to add a time-dependent pa- 
rameter a(t) in front of the fragmentation parameter (3{x). Then the alternance incubation-sonication 
corresponds to a rectangular function a(t) which is equal to 1 during the incubation time (since the 
sample is left to rest), and a max during the sonication pulse (where a max represents the maximal power 
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of the sonicator). We obtain the model 



d_ 

at' 



o 



u(x, t) = -V Q — (t(x)u{x, t)) - a(t)p(x)u(x, t) + 2a(t) / P{y)n{x, y) u{y, t) dy, (65) 



where u(x,t) still denotes the quantity of polymers of size x at time t. 

With this model, the problem of PMCA improvement becomes an optimization mathematical prob- 
lem: find a control a(t) which maximizes the quantity j xu(T, x) dx (total mass of pathogenic proteins) 
at a given final time T. The answer to this problem is difficult. First we can wonder if the rectangu- 
lar strategie used in the PMCA (alternance incubation-sonication) is the best one or if there exists 
a constant control which is better. This last question leads naturally to the problem of the fitness 
optimization for Equation ([65]) when a(t) = a is a time-independent parameter. Is a max the best 
constant to maximize A a ? Is there a compromise a op t G (1, a max ) to find? The answer depends on the 
coefficients r and f3 as indicated by the different theorems of this paper. More precisely, Theorem [2] 
ensures that the situation when an optimum a op t exists between 1 and a max can happen, and an 
example is given below. 



Example Let consider the same coefficients as in Figure 2(b) and suppose that the sonicator can 
multiply by 4 the fragmentation at its maximal power. Then in our model a max = 4 and we can see on 
Figure d] that the best strategy to maximize the fitness with a constant coefficient is not the maximal 
power but an intermediate a op t between 1 and a max - 




2 a op t 2 5 



> a 



Figure 4: The fitness is plotted as a function of a for the coefficients of Figure 2(b) There is 
a sonication value a op t in the interior of the window [1, a maa; ] which maximizes this fitness. 



In [7], the time dependent optimization problem is investigated on a discrete model. The study 
highlights the link between the optimal control a(t) and the constant a op t which optimizes the fitness 
of the system. 
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3.4 Therapeutic optimization for a cell population 

In the case when Problem (JT[) models the evolution of a size (or protein, or label, or parasites- 
structured cell population, r represents the growth rate of the cells and f3 their division rate. It is 
of deep interest to know how a change on these rates can affect the Malthus parameter of the total 
population, see for instance [HI [10]. It is possible to act on the growth rate by changing the nutrient 
properties - the richer the environment, the faster the growth rate of the cells. We can model such an 
influence by Equation ([3]), and the question is then: how to make A a as large (if we want to speed up 
the population growth, for instance for tissue regeneration) or as small (in the case of cancerous cells) 
as possible ? 

Plausible assumptions (see [H] for instance for the case of a size-structured population of E. Coli) for 
the growth of individual cells is that it is exponential up to a certain threshold, meaning that t{x) = tx 
in a neighbourhood of zero, and tending to a constant (or possibly vanishing) around infinity, meaning 
that the cells reach some maximal size or protein-content, leading to r —> < +oo. 

x— >oc 

Concerning the division rate f3, it is most generally vanishing around zero, either of the form f3{x) ~ 
fix 1 with 7 > 1 or with support [b, oo] with b > 0. It presents a maximum, and then decreases for large 
x - probably vanishing. Note that for r as for /3, very little is known about their precise behaviour for 
large sizes x, since such values are very rarely reached by cells in the real world. 

These assumptions allow us to apply our results. Theorem Q] and Corollary [1] lead to vanishing 
Malthus parameter A Q either for a — > or for a — > +oo. It means that against cancer, stressing the 
cells by diminishing nutrients can reveal efficient - this is very intuitive and it is known and used for 
tumors (by preventing them from vascularization for instance). What is less intuitive is that forcing 
them to grow too rapidly in size could also reveal an efficient strategy, as soon as it is established 
that the division rate decreases for large sizes (this last point could be studied by inverse problem 
techniques, see [331 E3 E] ) ■ It rejoins the same ideas as for Prions, as said in Section 1331 

On the opposite, in order to optimize tissue regeneration for instance, these results tend to prove 
that there exists an optimal value for a such that the Malthus parameter is maximum. This value 
can be established numerically (see Section 13.11 and [20J ) as soon as the shape of the division rate is 
known, for instance by the use of the previously-quoted inverse problem techniques. 

Conclusion 

The first motivation of our research was to investigate the dependency of the dominant eigenvalue of 
Problem (pQ) upon the coefficients /3 and r, since a first wrong intuition, based on simple cases, was that 
it should be monotonic (see [3(8]). By the use of a self-similar change of variables, we have explored 
the asymptotic behaviour of the first eigenvalue when fragmentation dominates the transport term 
or vice versa. This lead us to counter-examples where the eigenvalue depends on the coefficients in a 
non-monotonic way ; moreover, these counter-examples are far from being exotic and seem perfectly 
plausible in many applications, as shown in Section [3j A still open problem is thus to find what would 
be necessary and sufficient assumptions on r and /3, or still better on the ratio so that A a or A a 
be indeed monotonic with respect to a or a. 

Concerning our assumptions, a first sight at the statement of Theorem [1] gives the feeling that only 
the behaviour of the fragmentation rate f3 plays a role in the asymptotic behaviour of X a , and only 
the ratio - in the one of A a . This seems puzzling and counter-intuitive. In reality, things are not 
that simple: to ensure well-posedness of the eigenvalue problems ([3|) and ([!]), Assumptions (|26p and 
(I27|) strongly link r with /3, so that a dependency on /3 hides a dependency on r and vice versa. 
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Moreover, the mathematical techniques used here (moment estimates, multiplication by polynomial 
weights) force us to restrict ourselves to the space V of functions of polynomial growth or decay. The 
questions of how relaxing these (already almost optimal, as shown in |16j ) assumptions and how, if 
possible, express them in terms of pure comparison between r, k and f3 like in Assumptions (|26[) and 
(|27p are still open. 
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